Anomalous fluctuations of currents in Sinai-type random chains with strongly 

correlated disorder 



(N 

o 

(N 
00 



Gleb Oshanin/'Q Alberto Rosso, 2 -0 and Gregory Schehr 2 '0 

1 Laboratoire de Physique Theorique de la Matiere Condensee (UMR CNRS 7600), 
Universite Pierre et Mane Curie (Pans 6) - 4 Place Jussieu, 75252 Paris, France 
2 Laboratoire Physique Theorique et Modeles Statistiques (UMR CNRS 8626), Universite de Paris-Sud, Or say Cedex, France 

(Dated: September 19, 2012) 

We study properties of a random walk in a generalized Sinai model (SM), in which a quenched 
random potential is a trajectory of a fractional Brownian motion with arbitrary Hurst parameter 
H, < H < 1, so that the random force field displays strong spatial correlations. In this case, 
the disorder-average mean-square displacement (MSD) grows in proportion to log 2//H (n), n being 
time. We prove that moments of arbitrary order k of the steady-state current Jl through a finite 
segment of length L of such a chain decay as L~ (1 ~ H \ independently of k, which suggests that 
despite a logarithmic confinement the average current is much higher than its Fickian counterpart 
in homogeneous systems. Our results reveal a paradoxical behavior such that, for fixed n and 
L, the MSD decreases when one varies H from to 1, while the average current increases. This 
counter-intuitive behavior is explained via an analysis of representative realizations of disorder. 
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Since the pioneering works 1-3], random walks in ran- 
dom environments have attracted a great deal of atten- 
tion. In part, due to a general interest in dynamics in 
disordered systems, but also because such random walks 
have found many physical applications. Among them we 
mention dynamics of the boundary separating the coil 
and helix phases in a random heteropolymer [J, Q , dy- 
namics in a random- field Ising model [y, motion of 
dislocations in disordered crystals |, mechanical unzip- 
ping of DN A [9|1 , translocation of biomolecules through 
nanopores [10( and the dynamics of molecular motors 



111 ]. Some functionals arising in this model, like proba- 
bility currents in finite samples, show up in mathematical 
finance [H|,[l3|. Other examples can be found in [14- 16|. 

In the discrete formulation, one considers a random 
walk on a discrete lattice evolving in a discrete time. 
At each time step the walker jumps from site X to site 
X ± 1 with site-dependent probabilities px to X + 1 and 
q x = 1 - Px to X - 1 given by p x = \{l + e ■ s x ), 
where the amplitude < e < 1 measures the strength 
of the disorder and sx are independent and identically 
distributed (i.i.d.) random variables (r.v.). One often 
assumes binomial r.v., i.e., sx = ±1 with probability p 
and 1 — p. The probability P n (X) of finding the walker 
at site X after n time steps satisfies 

Pn(X) = q X+ l Pn-l(X + 1) + px-1 Pn-l(X - 1) . (1) 

In case of no global bias (p = 1/2), i.e. for the so- 
called Sinai model (SM), a remarkable result 0] is that 
for a given environment {px} the squared displacement 



Xf. 



miipx}) ln 4 (n) 



(2) 



as n — > oo with probability almost 1, where m({px}) is a 
function of the environment only [13] • Another intriguing 
feature of the SM concerns transport properties. It was 



revealed by analysing the probability current Jl through 
a finite Sinai chain of len gth L, that the disorder-average 
current decays as l/VL [18h21| . Curiously enough, de- 
spite a logarithmic confinement @, the disorder- average 
current appears to be anomalously high, so that such 
disordered chains offer on average less resistance with re- 
spect to transport of particles than homogeneous chains 
(all px = 1/2) for which one finds Fick's law Jl ~ 1/L. 

It is well-known that random walk in uncorrelated ran- 
dom environment {px} can be considered as a random 
walk in presence of a random potential Vl which repre- 
sents itself a random walk in space. Indeed, on scales of 
order L a random walk "explores" the potential 
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which is, clearly, a trajectory of a random walk with step 
length a = ki((l + e)/(l — e)). The case of standard Sinai 
diffusion where the r.v. sx are independent is now very 
well understood. On the contrary, there hardly exist an- 
alytical results in the case where these r.v. are strongly 
correlated. Such correlations are important, for instance, 
to understand the dynamics of the helix-coil boundary 
in random heteropolymcrs, where the location of chem- 
ical units usually exhibits strong correlations j22[. Such 
correlations are also currently studied in mathematical 
finance, as a necessary improvement over the standard 
Black-Sholes-Merton (BSM) model [H]. Any exact re- 
sult for such situations would thus be welcome and here 
we make a step further in this direction. 

In this Letter, we study properties of random walks in 
random environments in which the transition probabil- 
ities {px} are strongly correlated so that the potential 
Vl in j3]) is a fractional Brownian motion (fBm): Vl is 



2 



Gaussian, with Vl= 



and moments 



E{V L } = 0, E{(V l -V l ,) 2 }=<j 2 \L-L 



l\2H 



(4) 



where E{. . .} here and henceforth denotes averaging over 
realizations of Vl and < H < 1. The case H = 1/2 cor- 
responds to the original SM. For H < 1/2 the potential 
is subdiffusive while for H > 1/2 it is supcrdiffusivc. 

The mean squared displacement E{X^} in a correlated 
random environment can be estimated as follows. As- 
suming Arrhenius' law for the activated dynamics 
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the time required for a particle to diffuse in a disor- 
dered potential Vl over a scale of order of L, is of order 
of tx ~ e l , where V£ is a typical energy barrier. For 
Vl in (@|, V£ ~ crL H for typical realizations of disorder, 
so that one finds that for sufficiently large times, 

E{^}~a- 2 /"ln 2 / ff (n). (5) 

Our focal interest here is the understanding of the be- 
havior of the disorder-average current Jl through a fi- 
nite sample (of length L) of such a disordered chain, of 
its moments of arbitrary order, and eventually, of the full 
probability density function (pdf) of Jl- We proceed to 
show that, while its typical value is exponentially small 
J typ ~ exp (—L H ), its moments decay algebraically 



A k L-° , L> 1, 



(6) 



where 9 = 1 — H is the persistence exponent of the fBm 
2J, |25| and A k is a L-independent constant. We remind 



that the persistence exponent associated to a stochastic 
process characterizes the algebraic decay of its survival 
probability S(t) ~ t~ e " 
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The result in © is rather astonishing: it states first 
that the moments of the current of an arbitrary order 
decay in the same fashion. Second, one observes that for 
arbitrary H, < H < 1, the disorder-average current in 
such random chains is larger than the Fickian current in 
homogeneous systems. Finally, on comparing © and © 
for fixed n and L sufficiently large, and varying H, one 
concludes that E{X 2 } increases when H goes from 1 to 
0, while the disorder-average current decreases, which is 
an absolutely counter-intuitive and surprising behavior. 

In what follows we prove ([6]) and explain this aston- 
ishing behavior using three different, complementary ap- 
proaches: (i) a rigorous one, based on exact bounds, for 
the discrete random walk in a fBm potential, (ii) scaling 
arguments for the continuous version, both in space and 
time, of this diffusion (Langevin dynamics), which also 
allows to investigate the pdf of Jl (not only the moments) 
and (iii) thorough numerical simulations. We argue that 
(|6]) also holds for any disordered potential Vl in (J3]) which 
is the trajectory of a stochastic process with persistence 
exponent 9. As shown below, the occurrence of 9 in <j6j) 
can be physically understood as the moments of Jl are 
dominated by the configurations of Vl which drift to — oo 



without re-crossing the origin (see also [28jj). Such con- 
figurations occur with a probability ~ L~ e (L playing 
the role of time), yielding the L-dependence in ©. 




FIG. 1. (Color online) The disorder-average current E{Jl} 
(squares) and fj,k (L) with k — 2 (circles) and k — 3 (triangles) 
as a function of L for the fBm with H = 0.75, H = 0.4 and 
H = 0.25 (from top to bottom). The solid line is the predicted 
behavior ~ L~ e (|6]) with 6 = 1 — H. The temperature is 
T = 0.25 and averaging is performed over 10 5 samples. We 
used arbitrary units [a.u.] because the data for each value of 
H have been vertically shifted for clarity. 

We first consider the discrete version of SM, take a 
finite segment of length L and impose fixed concentra- 
tions of particles at the endpoints, P n (X = 0) = Pq and 
P n (X = L) = Pl- We then solve |T]) for a fixed envi- 
ronment {px}, and find that the steady-state current of 
particles through the chain Jl is given by 18, lj| 



Jl = 



DqPq _ DqPl 
tl r? 



(7) 



where Dq = 1/2 is the diffusion coefficient of a homoge- 
neous chain, tl is the so-called Kesten variable [2II : 



1 1 Pi 

T L = H 



P1P2 



qi qi q_i 



P1P2 ■ ■ -Pl-i 
qiq 2 ■ ■ ■ qL-i 



(8) 



and t£ is obtained from (JSJ) by replacements p k — > qL-k 
and q k — > PL-k- Thinking of L as "time", one notices 
that tl and r£ are time-averaged discretized geometric 
fractional Brownian motions (they can be thought of as 
the " prices" of Asian options within the framework of the 
fractional BSM model |13|). Note that in absence of a 



global bias E{1/tl} = E{1/t£}, and hence, without any 
loss of generality we set Pl = in what follows. Thus, 
combining ((H [3 [8]) and setting Pq = 1 yields 



Jl 




L-l 



X)exp(Vz) 



(9) 



1=1 



For typical realizations of {px}, the size of \Vi\ is 0(l H ) 
so that the typical current Jt yp is Jt yp ~ exp (~L H ) . We 
now turn to the disorder-average current and establish 
upper and lower bounds on fik{L) ©■ 
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To obtain an upper bound, consider a given realization 
of the sequence V\, V2, ■ ■ ■ , Vl—i and denote the maximal 
value among them as V max = maxo<i<L-i^- From © 
one has t l = (1 + J2i=i v l) > ex P (Knax), so that 



J k L < (1/2) K exp (-kV m 







(10) 



Hence, averaging in (| 10[) amounts to calculating the char- 
acteristic function of the distribution of Kn ax - Since 
exp (— fcVmax) vanishes in the limit L — > 00 (recall that 
Vmax ~ L H ) the average value of exp (— kV ma , x ) is dom- 
inated by configurations of the potential with V max — > 
0. The asymptotic behavior of the pdf Pl for 
fixed Vmax and large L is known from 2J, |25[ , yielding 
lnP L {V max ) = BlnL- 1 + 0{1), where 9 = 1 - H is the 
persistence exponent |3p| . Therefore, for L — > 00 one has 



Mfc (L) < B fe L 



(11) 



where B k is an L-independent constant. 

To determine a lower bound on Lik(L) we follow (l8l . 
19L 31 1 and make the following observation: averaging the 
expression in ([9]) amounts to taking average over the set 
ft of all possible trajectories {V/}i<;<l- Since tl > 0, a 
lower bound on ilk [L) can be straightforwardly obtained 
if one averages instead over some finite subset fl' C of 
trajectories with some prescribed properties: 



Ltk(L) > En> {J k L } 



(12) 



We choose ft' as comprising all possible trajectories 
{Vi}q<i<l which, starting at the origin at I = 0, never 
cross the deterministic curve Yi = Yq — odn(l + I) with 
Yq > and a > 1. For any such trajectory 



L-1 



L-1 



TL = l + ]Tcxp(^) < J2 

exp (Yi) < exp(Y Q )((a) , 
(13) 



z=i ;=o 
where C( a ) is the zeta- function, which implies that 
Hk(L) > E n - { Jt) > (cxp(yo)C(a)/2)- fe IEo' {1} . (14) 



One recognizes that Eq< {1} is, by definition, the sur- 
vival probability, Sl up to time L, for a fBm, starting 
at the origin in presence of a "moving trap" starting at 
Yo and evolving via Yi = Yo — aln(l + I). For standard 
Brownian motion (H — 1/2) in presence of a trap which 
moves as — l z , the leading largc-L behavior of the sur- 
vival probability Sl is exactly the same as in the case 
of an immobile trap, provided that z < 1/2 33]. It is 
thus physically plausible to suppose that the same be- 
havior holds for a more general Gaussian process such as 
a fBm. That is, one expects that for any H > the lead- 
ing large-L behavior of Eq> {1} will be exactly the same 
for an immobile trap and a logarithmically moving trap, 
i.e., that S L = E n - {1} - Y® /H /L e as L -> 00 0, El , 



where 9 = 1 — H (and this can be shown rigorously, up 
to logarithmic corrections pj|). Consequently, we find 



MfcW > D k L H - 



(15) 



where D k is independent of L. Note that the bounds 
in (fTUf and (fT5|) show the same L-dependence and thus 
yields the exact result announced in ([5]). 

We now turn to a continuous-time and -space descrip- 
tion of the dynamics in a disordered fBm potential. The 
position of the particle at time t is defined by x(t) € [0, L\ 
which obeys a Langevin equation 



X = -V'(x) +7)(t) 



(16) 



where V'(x) is a quenched random force [related to the 
hopping rates px in {J)] such that V(x) is a fBm with 
Hurst exponent H ((4]) and 77(f) is a Gaussian thermal 
noise of zero mean and covariance (r)(t)r](t')) = 2TS(t — 
f'). The current Jl can be obtained from the Fokker- 
Planck equation corresponding to (jTHJ) and reads 




exp[V(x)/T] dx 



(17) 



[see (0 © with D = T dS) and P = 1]. We remark 
here that the functional on the right-hand-side of (flTl) 
has been studied within a different context in [iil l32| . 
and it has been shown that its first moment decreases as 
L~ 6 with 9 = 1 — H. We focus next on the moments of 
higher order and on the pdf of the functional in (jT7j) . 

Instead of Jl/T, which can be viewed as the in- 
verse of the partition function of SM, we focus on 
the free energy F = T log( Jl/T) and study its pdf 
Ht=o(F). Let us first discuss the case T = where 
= mino<;c<i, V(x). Recalling that V(0) = 
(such that i? m i n < 0), the cumulative distribution 
qL(E) = Pv(E m m > —E), with E > 0, coincides with the 
survival probability of V(x) up to 'time' L starting in E 
at initial 'time' in presence of an absorbing boundary at 
V = 0. For fBm, and more generally for any sclf-affinc 
process, qL(E) takes the scaling form qi, (E) = Q(E/L H ). 
For L ^s> E x t H , gj,{E) behaves algebraically with a per- 
sistence exponent 6 Hfjj], q L (E) ~ E e / H / L e (9 = 1 - H 



for fBm [2J,l25|). In the opposite limit where L < E x / H , 



qL{E) decays as a Gaussian [for Gaussian V(a;)]. There- 
fore, one has for n T=0 (/) = d E qL{E)\ E= _ F 

, F > 

Il T=0 (F) = { L-^F] 9 / 11 - 1 , -L H « F < (18) 
cxp(~F 2 /2L 2H ) , F < -L H . 

What happens at finite T where F = E — TS is now 
the balance between the energy E and the entropy S ? 
One expects the particle to be localized close to the mini- 
mum E m i n , which is of order 0(L H ), while the maximum 
entropy ~ O(lnL). Hence when L 1 the main contri- 
bution to F comes from E m i n so that for a given sample 
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FIG. 2. (Color online) Pdf of the free energy ~F/L H for 
different system sizes: L — 64, 256, 1024 from bottom to top. 
The black curve represents the distribution of Emi n /L H for 
L — 4096. Inset: Pdf of £7 m i n — F for different system sizes: 
L = 64, 256, 1024 from left to right. Histograms are computed 
using 10° samples and setting T = 1. 



at finite T, F will be very close to E m - m . This is corrob- 
orated by our numerical simulations (sec Fig. [2]) ■ 

We now come back to the current distribution. Very 
small currents, Jl <C Jt yp ~ exp(— L H ), correspond to 
F < -L H in (JT5J) and one obtains that P{Jl) is log- 
normal, InP(Jk) cx — ln 2 (Ji/T). Within the opposite 
limit, Jl S> Jtypj one finds from (p~8|) that 



P{Jl) 



pog(J £ /T)] 



e/H-i 



(19) 



This algebraic behavior holds up to a large cut-off value 



Jmax- At T = we have a sharp cut-off at J max = 1 (flgf . 
at finite T, P(Jl) has a fast decay which depends on the 
fluctuations of V(x) at a short length scale close to the 
origin x = 0. For L > 1, the moments of the current 
are dominated by the regime where Jt yp <C Jl < Jmax ~ 
0(1) (H9|) such that one gets fx k (L) ~ l/L 9 , ©. 

To perform numerical simulations we consider a dis- 
crete random potential Vfc, fc = 0, 1,...L — 1, which 
displays fBm correlations (J4|, we set tr 2 = 1. We 
use a powerful algorithm, as explained in J^H l3o| . 
which allows to generate very long samples of fBm 
paths. For each sample, we compute the current, Jl = 
T lEkZo cM~V k /T)]-\ the free energy F = T\n( J L /T) 
and the ground state energy E m i B = min^ 14. In Fig. [1] 
we show a plot of the first three moments as a function 
of L for different values of H. These plots show a very 
good agreement with our analytical predictions in ([6]). In 
Fig. [2j we show that the pdf of the rescaled free energy, 
F/L H , at finite temperature T converges to the pdf of 
the rescaled ground state energy E min /L H . The reason 
for this is that, for each sample, the difference between 
F and E m i n grows very slowly with L, probably logarith- 



mically (inset of Fig. [2]). In the rescaled variables, this 
difference vanishes when L 



CO. 




FIG. 3. (Color online) The pdf in Eq. (21) for o 2 L 2H = 1/2 
(blue), a 2 L 2H = 2 (red) and a 2 L 2H = 5 (green). 

We close this paper with an observation that such 
Sinai-type chains show an interesting transition to a 
diode-like behavior when a (at fixed L) or L (at fixed a) 
exceeds some critical values. To illustrate such a transi- 
tion, consider a finite chain in which at site X — we 
maintain a fixed concentration Pq = 1 of one type of 
"white" particles and place a sink for them at X = L. 
On the other hand, at X = L we introduce a source 
which maintains concentration 1 of "black" particles and 
place a sink for them at site X = 0. The particles are 
noninteracting and in the steady-state, for fixed {px} we 
have two counter-currents of white (J™) and black ( j£) 
particles. Obviously one has E{(J L ") k } = E {(j£) fc } for 
any k > 0. Consider next the random variable 



W = J L U /(J L U + Jl)=T* L /(T* L +T L ), 



(20) 



which probes the likelihood of an event that for a given 
realization of disorder J^' = J h L . The pdf of u> can be 
calculated exactly to give 



P( W ) 



1 



In 2 



2ttw(1 - uj)crL H 



exp 



2<r 2 L 2H 



(21) 



A remarkable feature of (|2T|) is depicted in Fig. El 
Namely, P(w) changes modality when o x l H L crosses a 
critical value a¥ H L c = 2 1 / 2H . For short chains (small 
disorder) P(ui) is unimodal and centered at oj = 1/2: any 
given sample is transmitting particles in both directions 
equally well and, most probably, J^ = J b L . On contrary, 
for sufficiently strong disorder and/or long chain, P(co) 
becomes bimodal with a local minimum at uj = 1/2 and 
two maxima close to and 1. This means that a given 
sample is most likely permeable only in one direction. 
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